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Abstract 

Using a dataset corresponding to an integrated luminosity of 3.0 fb _1 collected in pp 
collisions at centre-of-mass energies of 7 and 8 TeV, the B® —> (pep branching fraction 
is measured to be 

B(B° S —> (pep) = (1.84 ± 0.05 (stat) ± 0.07 (syst) ± 0.11 (f s /fd) ± 0.12 (norm) )xl0~ 5 , 

where f s /fd represents the ratio of the B® to B° production cross-sections, and the 
B° —> <f>K*( 892)° decay mode is used for normalization. This is the most precise 
measurement of this branching fraction to date, representing a factor five reduction 
in the statistical uncertainty compared with the previous best measurement. A 
search for the decay B° —> epep is also made. No signal is observed, and an upper 
limit on the branching fraction is set as 

B(B° -> (pep) < 2.8 x 1(T 8 

at 90% confidence level. This is a factor of seven improvement compared to the 
previous best limit. 
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1 Introduction 


In the Standard Model, the flavour-changing neutral current decay B ° — y 00 proceeds via 
a b — > sss penguin amplitude. The decay was first observed by the CDF experiment at the 
Tevatron (!]. Subsequently, it has been studied by the CDF and LHCb collaborations, who 
searched for CP-violating asymmetries in the decay time and angular distributions of this 
mode |2--j5j. These studies provide a probe for possible new physics contributions entering 
into the penguin loop and B® — B° s mixing diagrams (6|. Furthermore, as the B° s — > 00 mode 
will be used as normalization for studies of other charmless B® meson decays, it is important 
to have a precise determination of its branching fraction. The CDF collaboration measured 
this relative to the decay B® —>• J/00 |2|. Using the current value of the B® —» J/00 
branching fraction 0 , the CDF result gives B(B° S —» 00) = (1.91 ± 0.26 ± 0.16) x 10 5 , 
where the first uncertainty is from the measured ratio to B° s —)■ J/-00, and the second is 
due to the knowledge of the B° s —>■ J/00 branching fraction. Various predictions from 
theories based on QCD factorization exist for the B® —> 00 branching fraction (18-10 


These suffer from uncertainties related to weak annihilation diagrams. These uncertainties 
are controlled using experimental information from decays such as B° —> 0R+892) 0 . 
Several recent predictions are summarized in Table [l] The central values are in the range 
(1.5 - 2.0) x 10+ 

In this paper the B° s —00 branching fraction (the use of charge-conjugate modes is 
implied throughout) is measured using the full LHCb Run 1 dataset, comprising data 
corresponding to an integrated luminosity of 1.0 fb” 1 collected in pp collisions at a centre- 
of-mass energy of 7 TeV, and 2.0fb -1 collected at 8 TeV. The decay B° —> 0/1 *(892)°, 
which has a similar topology, is used for normalization. The 0 and /i*(892)° mesons 


are reconstructed in the K + K~ and K + n 
search for the yet unobserved decay B° 
Standard Model by the OZI rule 12 


final states, respectively. In addition, a 
—> 00 is made. This decay is suppressed in the 
with an expected branching fraction in the range 


(0.1 — 3.0) x 10~ 8 [8,10,13,14]. However, the branching fraction can be enhanced, up to 
the 10 - ' level, in models such as supersymmetry with R-parity violation 14 . The current 
best limit for this mode is from the BaBar collaboration |15|, B(B° —> 00) < 2.0 x 10“' 
at 90% confidence level. 


Table 1: Predictions for the B® —> 00 branching fraction. The first and second uncertainties of 
Refs. [8,9 reflect the knowledge of CKM parameters and power corrections, respectively. 


B(B<> #) (10- 5 ) 

Approach 

Reference 

1.95 ± 0.10«;S 

i ^7+0.26+1.13 
± ‘ u * -0.21-0.88 

1 e ; e;+ 2 - 24 

1 «7+ 0 - 89 

1 - D ‘ -0.71 

QCD factorization 
QCD factorization 
QCD factorization 
pQCD 

I 

8; 

9 

10 
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2 Detector and software 


The LHCb detector |16fl7| is a single-arm forward spectrometer covering the pseudorapidity 
range 2 < 77 < 5, designed for the study of particles containing b or c quarks. The detector 
includes a high-precision tracking system consisting of a silicon-strip vertex detector 
surrounding the pp interaction region 118], a large-area silicon-strip detector located 
upstream of a dipole magnet with a bending power of about 4 Tm, and three stations of 
silicon-strip detectors and straw drift tubes 19 placed downstream of the magnet. The 
tracking system provides a measurement of momentum, p, with a relative uncertainty that 
varies from 0.5% at low momentum to 1.0% at 200GeV/c. The minimum distance of a 
track to a primary pp interaction vertex (PV), the impact parameter, is measured with a 
resolution of (15 + 29/px) M- m , where px is the component of the momentum transverse to 
the beam, in GeV/c. 

Different types of charged hadrons are distinguished using information from two 
ring-imaging Cherenkov detectors 1201. Photon, electron and hadron candidates are 
identified by a calorimeter system consisting of scintillating-pad and preshower detectors, 
an electromagnetic calorimeter and a hadronic calorimeter. Muons are identified by a 
system composed of alternating layers of iron and multiwire proportional chambers 


21 


The trigger 22 consists of a hardware stage, based on information from the calorimeter 
and muon systems, followed by a software stage, which applies a full event reconstruction. 
The software trigger applied in this analysis requires a two-, three- or four-track secondary 
vertex with a significant displacement from any PV. At least one charged particle must- 
have a transverse momentum px > 1-7 GeV/c and be inconsistent with originating from 
a PV. A multivariate algorithm |]23] is used for the identification of secondary vertices 
consistent with the decay of a b hadron. 

In the simulation, pp collisions are generated using Pythia 
configuration 25 . Decays of hadronic particles are described by EvtGen 
final-state radiation is generated using PHOTOS [27]. The interaction of the generated 


with a specific LHCb 
in which 


26 


particles with the detector, and its response, are implemented using the Geant4 toolkit 
as described in Ref. 
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3 Signal selection 

The selection of candidates takes place in two stages. First, a selection using loose criteria 
is performed that reduces background whilst retaining high signal efficiency. Following 
this, a multivariate method is used to further improve the signal significance. 

The selection starts from charged particle tracks that traverse the entire spectrometer. 
Selected particles are required to have px > 500 MeV/c. Fake tracks created by the 
reconstruction due to random combinations of hits in the detector are suppressed using 
a requirement on a neural network trained to discriminate between these and genuine 
tracks associated to particles. Combinatorial background from hadrons originating at the 
primary vertex is suppressed by requiring that all tracks are significantly displaced from 
any primary vertex. Kaon and pion candidates are selected using the information provided 
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by the ring-imaging Cherenkov detectors. This is combined with kinematic information 
using a neural network to provide an effective probability that a particle is a kaon (V K ) 
or pion ( V 7r ). To select kaon candidates it is required that V K (1 — V n ) > 0.025. The pion 
candidate in the B° —> 0/1 *(892)° decay mode is required to have V n > 0.2 and V K < 0.2. 

The selected charged particles are combined to form 0 and K* meson candidates. The 
invariant mass of the K + K~ (/i + 7r _ ) pair is required to be within 15MeV/c 2 (150MeV/c 2 ) 
of the known mass of the 0 (/l*(892)°) meson |7|. In addition, the p t of the 0 and K* 
mesons must be greater than 1 GeV/c. 

Candidates for the decay B° s —1 (jxj) are formed by combining pairs of 0 mesons. A 
fit is made requiring all four final-state particles to originate from a common vertex, 
and the direction vector between the primary and secondary vertices is required to be 
consistent with the direction of the momentum vector of the B° s meson candidate. Further 
requirements are then applied to remove background from specific 6-hadron decays that 
peak close to the B® mass. To reject background from B° —i 0/7* (892)° decays, the kaon 
with the lowest value of V K is considered to be a pion, and the K + n~ and Ii + K~K + n~ 
invariant masses are calculated. Candidates with m(K + n~) within 50MeV/c 2 of the known 
K* (892)° mass and m(K + K~K + TT~) within 30MeV/c 2 of the B° mass |7| are rejected. 
Similarly, to remove decays via open charm mesons, the K + K~ir + mass is calculated. If 
m(K + K~ 7r + ) is within 22.5MeV/c 2 of the D + or Df mass |7|, the candidate is rejected. 
These vetoes are found to retain 91% of simulated B ° —> 00 decays. 

Candidates for the decay B° —> 0/1* (892)° are formed from combinations of 0 and K* 
mesons. Identical vertex and pointing requirements as for the B° s —> 00 decay mode are 
applied. To reject background from B° s —> 00, the mass of the K + pair is calculated 
assuming that both hadrons are kaons. Candidates with m(K + K~) within 15MeV/c 2 
of the 0 mass and m(K + K~K + K~) within 30MeV/c 2 of the B° s mass |7| are rejected. 
Background from open charm decays is suppressed in a manner similar to that used for the 
B° s —> 00 candidates. These vetoes are found to retain 97% of simulated B° — * 0/T*(892)° 
decays. 

The combinatorial background is further suppressed using a Boosted Decision Tree 
method (BDT) |30,31 . The BDT is trained to identify four-body hadronic 6-hadron decays 
with high efficiency using independent data samples of such decays. It uses information on 
the displacement of the 6-hadron candidate from the primary vertex, kinematic information 
and track isolation criteria. Although the same BDT is used for the B® —> 00 branching 
fraction measurement and the search for B° —y 00, the method used to optimize the cut on 
the BDT output is different. For the branching fraction measurement, the cut optimization 
is based on the normalisation mode B° —>■ <f>K*( 892)°. The figure of merit used is 

So x £3 


\/Sq x e s + fV bg 


where So is the signal yield of B° —y <pK*( 892)° candidates in data before any BDT cut is 
applied, £s is the efficiency of the BDT cut on simulated B° —> 0/1 *(892)° decays, and 
A0 g is the number of background candidates surviving the BDT cut in a suitable upper 
sideband of the 0/i*(892)° candidate mass distribution, scaled to the width of the B° 
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signal window. Maximizing this figure of merit results in a rather loose BDT requirement 
that retains 98% of signal events while rejecting more than 90% of the background. 

For the B° —> 00 search, the figure of merit used is 


a / 2 + 


with a set to 3, corresponding to the signal significance required to claim evidence for a 
new decay mode 


32 . Here e' s is the efficiency of the BDT cut on simulated B° s —> 


decays, and is the number of background candidates surviving the BDT cut in an 
upper sideband of the 00 candidate mass distribution, scaled to the width of the B® signal 
window. Maximizing this figure of merit results in a tighter BDT requirement that retains 
87% of signal events. 


4 Fits to mass spectra 

The yields for the signal and normalization channels are determined from fits to the 
invariant mass distributions of the selected candidates. In the simulation, the B° s — > 00 
invariant mass distribution is well modelled by a probability density function (PDF) 
consisting of the sum of three Gaussian distributions with a common mean. In the fit- 
to the data, the relative fractions of the Gaussian components are fixed to the values 
obtained from the simulation, whilst the widths are allowed to vary by an overall resolution 
scale factor. The yield and common mean are also left free. After applying all selection 
requirements, the only remaining background is combinatorial, which is modelled by a 
constant. No component for B° —>- 00 decays is included in this fit. Figure [l] shows the 
resulting fit to data, which gives a signal yield of 2309 ± 49 candidates. 

The B° — y 0/1 *(892)° invariant mass distribution is modelled by a PDF consisting of 
the sum of a Crystal Ball function |33j and two Gaussian functions. As for the signal 
mode, the relative fractions of the components and the tail parameters are fixed in the fit 
to the data, whilst the widths are allowed to vary by an overall resolution scale factor. 
The yield and mean are also left free. A component is also included to account for the 
small contribution from the decay B ° —> 0/1 *(892)° [34j. The shape parameters for this 
component are shared with the B° component, while the relative position is fixed to the 
known mass difference between the B° and B® mesons |7|. Combinatorial background is 
modelled by an exponential function. 

Potential peaking backgrounds, from T° —>• 0p7r - (with the proton misidentified as a 
kaon) or A® — y <ppK~ (with the proton misidentified as a pion), are modelled using a single 
histogram PDF generated from simulated events. The relative yield of each decay mode is 
weighted according to the expectation from the simulation. The yield of this component 
is left to float- in the fit. Backgrounds from B® —> 00 and open charm decay modes are 
negligible after the vetoes described in Section [3] have been applied. Figure [2] shows the 
result of the fit of this model to the B° —> 0A'*(892)° dataset after all selection criteria 
are applied. The yield of B° candidates determined by the fit is 6680 ± 86. 
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Figure 1: The K + K~ K + K~ invariant mass distribution. The total fitted function as described 
in the text is shown by the (red) solid line, the B® —> (jxj) component by the (blue) long-dashed 
line, and the combinatorial background as the (purple) dotted line. 


5 Branching fraction for B° s -»■ U 


The branching fraction of B° s —y 
determined using 

B(B° S 00) 


relative to that of the B° —> 0A"*(892)° decay mode is 


N h ^1* ( 892)o B(K* (892)° -)■ K+tT 


B(B°->0AT*(892)0) N<pK , {892)0 e 


sel 


B{<t> K+K- 


fs/fd 


where the terms B are the branching fractions of the stated decay modes, N are the 
signal yields, £ sel are the selection efficiencies, and the fragmentation fraction ratio, 
fs/fd, is the ratio of the B® to B° production cross-sections. The selection efficiencies 
are determined from simulation, apart from those related to the particle identification, 
which are determined in data using large calibration samples of charged kaons and pions 
from D* + —y D°(—} K ~ti + )'k + decays 20 . The ratio of efficiencies is found to be 


e 


sel 


892)°/ 


r sel 


/ e M) = 0.795 ± 0.007, where the uncertainty is purely statistical. The value of 


fs/fd is taken from previous LHCb analyses as 0.259 ± 0.015 35 -37 


The signal yields are determined using the mass fits described in Section |4} These 
values are corrected for the fraction of candidates where one of the hadron pairs, K + K~ or 
K + 7T — , is produced in a non-resonant S-wave configuration, rather than as a 0 or Jl*(892)°. 
The S-wave fractions are taken from previous LHCb angular analyses of the B® —> (jxj) and 
B° —> <f)K* (892)° decay modes. For the B° s —» 00 decay mode, we use the measured value 
of 2.1 ± 1.6 % [H] as the S-wave fraction within the K + K~ invariant mass range used for 
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Figure 2: The K + K~K + ir~ invariant mass distribution. The total fitted function is shown by the 
(red) solid line, the B° — > 4>K* component by the (blue) short-dashed line, the B® — > 0A*(892)° 
component by the (blue) long-dashed line, the A® —> (j)pn~ and —> <f>pK~ contribution by the 
(green) dashed-dotted line, and the combinatoral background by the (purple) dotted line. 


this analysis. Similarly, for the B° —> 0A"*(892)° decay mode, we use a measured value 
for the S-wave fraction. The uncertainties on these fractions lead to 


To 38 


of 26.5 ± 1.8' 

a 3.1 % relative uncertainty on the ratio of branching fractions. This procedure assumes 
that the efficiencies for the P- and S-wave components are the same. In the simulation, a 
1.1 % difference is observed between these efficiencies, and this is assigned as an additional 
uncertainty. 

Various other uncertainties arise on the measurement of the ratio of branching fractions. 
The limited size of the available simulation samples leads to a relative uncertainty of 0.8 %. 
The influence of the assumed mass model is probed by performing the fit with different 
models for the signal and background components. This includes quantifying the effect 
of removing the peaking background component in the B° —> <j>K*( 892)° fit. The largest 
variation in the ratio of branching fractions seen in these studies is 0.6 %, which is assigned 
as a relative systematic uncertainty. 

The track reconstruction efficiency agrees between data and simulation at the level of 
2.0% 39 . This uncertainty largely cancels in the ratio of branching fractions. A residual 
relative uncertainty of 0.5% remains due to the fact that the pion in the B° —>■ 0/1 *(892)° 
decay mode is relatively soft. An additional relative uncertainty of 0.3% is assigned to 
account for the difference in the hadronic interaction probabilities for kaons and pions 
between data and simulation. A further uncertainty arises from the modelling of the 
hardware trigger in the simulation. This is estimated using a data-driven technique and 
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Table 2: Summary of the systematic uncertainties on the measurement of the ratio of branching 
fractions B(B® (jxj))/B(B 0 —> 4>K*). 


Source of systematic uncertainty 

Relative uncertainty (%) 

S-wave fraction 

3.1 

Relative efficiency between P and S-wave 

1.1 

Simulation sample size 

0.8 

Fit model 

0.6 

Tracking efficiency 

0.5 

Hadronic interactions 

0.3 

Hardware trigger 

1.1 

Particle identification efficiency 

0.3 

H(0 -)• K+K~) 

1.0 

Quadratic sum of the above 

3.8 

Fragmentation fraction ratio ( f s /fd ) 

5.8 


leads to a relative systematic uncertainty of 1.1% on the ratio of branching fractions. 
Variations in the procedure used to determine the relative particle identification efficiency 
lead to a relative uncertainty of 0.3 %. Possible systematic effects on the efficiency for 
Bg —> 00 due to the finite width difference in the B° s system 40 have been checked, and 
found to be negligible. The value of B{4> —>■ K + K ~) is taken from Ref. |7| and contributes 
a relative uncertainty of 1.0%. The value of B(K*( 892)° —y K + n~) is taken to be 2/3 
exactly. The systematic uncertainties are summarized in Table [2j Summing these in 
quadrature gives a relative uncertainty of 3.8% on the ratio of branching fractions. The 
knowledge of the fragmentation fraction ratio, f s /fd, gives a relative uncertainty of 5.8%, 
which is quoted separately. 

The ratio of branching fractions is found to be 


B(B° a -> 00) 
B(B° 0K*) 


1.84 ± 0.05 (stat) ± 0.07 (syst) ± 0.11 (f s /fd)- 


This is converted into an absolute branching fraction using B(B° —> 4>K*( 892)°) = 
(1.00 ±0.04 ±0.05) x 10 -5 , which is obtained by averaging the results in Refs. j4l] and [42 
assuming that the uncertainties due to the fragmentation fractions and S-waves are fully 
correlated between the two measurements. The resulting value for the absolute branching 
fraction is 


B(B° —> 00) = (1.84 ± 0.05 (stat) ± 0.07 (syst) ± 0.11 (f s /fd) ± 0.12 (norm)) x 10 5 . 
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6 Search for the decay B° -»• 4>(f> 


To search for the B° —> 00 decay mode, the tight BDT selection described in Section [ 3 ] is 
used. To fit for a putative B° —y 00 signal, the same signal model as for the B ° signal 
is used. The mean value of the signal mass is shifted relative to the B® mode by the 
known B®-B° mass splitting, and the resolution parameters are kept common between 
the two modes. The resulting fit is shown in Figure |3j The data are consistent with 
having no B° —> 00 contribution. The fitted B° signal has a yield of 5 ± 6 events, and the 
statistical significance is less than 2 standard deviations, hence an upper limit is placed on 
the branching fraction of the decay. 

To determine this limit, a modified frequentist approach, the CL S method, is used 
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The method provides CL s+ b, a measure of the compatibility of the observed distribution 
with the signal plus background hypothesis, CLb, a measure of the compatibility with 
the background only hypothesis, and CL S = CL s+b /CL b . The expected and observed CL S 
values as a function of the branching fraction are shown in Figure [4j This gives, at 90% 
confidence level, an upper limit of B(B° —» 00) < 2.8 x 10~ 8 . At 95% confidence level, 
the upper limit is found to be B(B° —> 00) < 3.4 x 10 -8 . 



Figure 3: The K + K~K + K~ invariant mass with the tight BDT selection applied. A fit to the 
total PDF as described in the text is shown as a (red) solid line, B® —> 00 as a (blue) long-dashed 
line, B° —> 00 as a (blue) short-dashed line, and the combinatorial background as a (purple) 
dotted line. 
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Figure 4: Results of the CL S scan as a function of the B° —> 00 branching fraction ( BF ). 
The observed CL S distribution is given by the (black) points and solid line, while the expected 
distribution is given by the (black) dashed line. The dark (green) and light (yellow) bands mark 
the la and 2a confidence regions on the expected CL S . The upper limits at 90% and 95% 
confidence level are where the observed CL S line intercepts the (red) solid and dashed horizontal 
lines, respectively. 


7 Summary 


The ratio of branching fractions B(B® —> 00)/i3(i? 0 —> 0/1* ) is determined to be 

= 1-84 ± °'° 5 ( Stat ^ ± °'° 7 ^ SySt ) ± °' n 

where the first uncertainty is statistical, the second is systematic, and the third is due 
to the ratio of fragmentation fractions. The absolute branching fraction for B G S —> 00 is 
determined to be 


B(B ° —y 00) = (1.84 ± 0.05 (stat) ± 0.07 (syst) ± 0.11 {f s /fd) ± 0.12 (norm)) x 10 5 . 


This is in agreement with, but more precise than, the measurement made by the CDF 
collaboration, B(B° S —> 00) = (1.91 ± 0.26 ± 0.16) x 10 -5 . It is also in agreement with 
theory predictions [8- 

A search for the decay B° —> 00 is also made. No significant signal is seen, and an 
upper limit of 
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B(B° -4 00) < 2.8 x 10" 8 


is set at 90% confidence level. This is more stringent than the previous limit of 
B(B° —>■ 00) < 2.0 x 10” 7 , set by BaBar 15 , and provides a strong constraint on possible 


contributions to this mode from physics beyond the Standard Model 14 
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